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Abstract 

The present work is essentially concerned with the development of statistical theory 
for the low temperature dislocation glide in concentrated solid solutions where atom- 
sized obstacles impede plastic flow. In connection with such a problem, we compute 
analytically the external force required to drag an elastic string along a discrete two- 
dimensional square lattice, where some obstacles have been randomly distributed. 
The corresponding numerical simulations allow us to demonstrate a remarkable 
agreement between simulations and theory for an obstacle density ranging from 1 
to 50 % and for lattices with different aspect ratios. The theory proves efficient on 
the condition that the obstacle-chain interaction remains sufficiently weak compared 
to the string stiffness. 
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1 Prom the solid solution strengthening theory 



The statistical theory for solid solution hardening (SSH) emerged from the 
seminal works of Sir N. Mott[T] and his near colleagues, F.R.N. Nabarro [2|l3] 
and J. Friedel The early analytical theory, perfected and extended by other 
contributors, as for instance R. Fleischer, R. Labusch and T. Suzuki [5|6|7|l8|9] 
applies to substitutional alloys where the solute atoms can be considered as 
immobile during the dislocation glide, by contrast to the cases where dislo- 
cations may drag along an atmosphere of fast diffusing impurities. In face 
centered cubic (fee) alloys, the critical resolved shear stress (CRSS) was then 
expected to increase in proportion to with Cs as the atomic concentration 
of solute atoms and r] as an exponent depending on the assumptions made on 
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the interaction between dislocations and foreign atoms: 77 = 1/2 in Friedel- 
Fleischer (FF) theory while 77 = 2/3 in Mott-Nabarro-Labusch (MNL) theory 
and ?7 = 1 in Friedel-Mott- Suzuki (FMS) [US]. Within analytical theory for 
SSH, the dislocation is thought of as a continuous elastic string impinged on 
a two-dimensional (2D) random static potential. The depinning transition in 
such a model is a typical issue of statistical physics, belonging to a broad 
class of problems concerned with extended interfaces motion in heterogeneous 
materials [T0][TT][TMM1[T5][TM71[TMM20^ . 

The recent developments of three-dimensional atomistic simulations (3D- AS) 
allowed to work on more realistic models for dislocations in solid Solutions |23f24ll25ll26f27ll28ll29f30l 
Though 3D-AS confirmed that a large part of the dislocation pinning hinges 
on the impurities situated in the crystal planes that bounds the dislocation 
glide plane [23], the simulations revealed also the complexity of the dislocation- 
obstacle interaction. In fee alloys, the geometry of the dislocation core, dis- 
sociated in two Shockley partials separated by a (111) stacking fault ribbon 
undermines the simple picture of an elastic line in interaction with a single 
type of obstacles, as stated in the basic version of SSH theory. Instead, the 
pinning forces differ according to partials and to the obstacle positions, i.e., 
above or below the glide plane [2SI2S]- 

On the other hand, the nanometric scale of the atomistic simulations, a strin- 
gent limit imposed by the computational load, hinders the direct extrapola- 
tion of simulation results to macroscopic samples. A multi-scale approach is 
thence required to link the atomistic studies to the realm of materials sciences. 
A manner to proceed consists in incorporating some of the atomic ingredients 
revealed in 3D-AS to a discrete version of the elastic string model that remains 
tractable even for large dimensions [31]. The discretization of the string model 
should allow to transfer of the atomic details. 



In the present paper, the impact of the discretization on the depinning transi- 
tion is analyzed thoroughly. The elastic string is replaced with a discrete spring 
chain, the nodes of which move on a 2D square lattice and interact with some 
pinning points randomly distributed on lattice sites. This very simple model 
allows us to devise an analytical theory which accounts for the discreteness of 
the obstacle distribution and thus opens a promising perspective to integrate 
more of the atomic details. In order to demonstrate the accuracy of the theory, 
we compute directly the critical external force within numerical simulations 
applied to the discrete string model. Theory and simulations agree remark- 
ably well for a broad range of model parameters, e.g., (i) the in-plane obstacle 
density c^, (ii) the lattice size in every direction of space, (iii) the maximum 
pinning force /m and (iv) the potential interaction cutoff w that characterizes 
the obstacles. The theoretical predictions are proved reliable on the condition 
that /m and w remain smaller than certain bounds varying with c^. For a dense 
distribution of weak pinning points, the critical configuration of the chain is 
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found to be a quasi-straight line parallel to the atomic rows. The depinning is 
then shown to occur at some vacant site clusters (VSC) which the typical size 
is explicitly related to the lattice dimensions in both directions of space. The 
external force required to drag along the spring chain over a finite distance 
reflects such a size dependence. Noteworthily the effective density exponent 
rj is also found to vary with lattice dimensions, in contrast with expectations 
drawn on standard SSH theory. 

Our report is organized as follows. In Sec. II, the spring chain model is intro- 
duced and the direct numerical computations are described. In Sec. Ill, the 
statistical theory is derived and compared to the numerical data. The results 
are resumed and commented in Sec. IV. 



2 The phenomenological spring chain model 

The model proposed hereafter belongs to the wide class of elastic interface 
models, extensively studied in statistical physics |16fl7fl8fl9f20ll21f22j . A one 
dimensional elastic string is discretized with a spatial step b, equivalent to 
the shortest interatomic distance in solids. Each node of the discrete chain 
is bound to its first neighbor by an harmonic spring of strength F. The two 
quantities, b and F are chosen to scale distances and forces, respectively. The 
spring chain nodes move along the column of a square lattice. The size of 
the lattice in the direction of the chain is denoted as Ly whereas the distance 
over which the chain is dragged is L^. The 2D random array of obstacles is 
constructed by selecting the occupied lattice sites, up to a number of obstacle 
equals to CsL^^Ly, where Cs is the obstacle density. Since the depinning process 
occurs when the chain nodes pass the force maximum, the interaction potential 
is expended as a polynomial function in the vicinity of such a maximum. 
Assuming that the interaction is attractive and that the potential is symmetric 
with respect to its minimum, we obtain a polynomial function of at least fourth 
order: 

V{x) = Vo{x'^/w'^ - 1)^ for |a;| < w 

V{x) =0 for |x| > w, (1) 

which corresponds to a force f{x) = — 4Vo(x^/tf^ — l)x/w'^, with a maximum 
value /m = 8| Vo|/(3\/3w), attained when x = ±w/y/3. The chain nodes inter- 
act solely with obstacles situated in the column along which they may glide. 
The polynomial force with a distance cutoff w is obviously very far from the 
dislocation-solute interaction, characterized by a decrease of Coulomb type. 
Hereby we describe only the local potential variation yielded when a solute 
atom visits a dislocation core. The parameter w fixes how the interaction de- 
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position of the chain / lattice X-axis 



Fig. 1. (color online) Strongest pinning configuration of the spring chain on a random array of obstacles 
(circles) for fm = 0.1, Ly = 2000, Lx = 500, a density Cs = 16 % and an interaction cutoff w = 1. 
Only obstacles close from the chain have been reported for clarity. Inset shows a magnification of obstacles 
(circles) and nodes (triangles) of a chain segment. X and Y axis have different scaling for convenience of 
the plot. 

creases in the vicinity of the force maximum. Both J'm and w can be extracted 
from atomistic data as those reported in [25||28] . 

The dimensionless over-damped Langevin dynamics for the chain node k is 
given by: 

Xk = [Xk+i + - 2xk] + T-Y: 4^ (^fe - Sk,) ^ {Xk - Sk,y _ 

~^ 

where is the position of the node /c, r is the external force and s^^i is the 
coordinate of the \th obstacle in the \.th row. For the weak pinning forces 
we are concerned with, the chain strain remains very small such that the 
anharmonic terms in the spring tension have been neglected. Properly scaled, 
the continuous version of the spring chain model served in the development of 
the SSH theory [T1I31I1I8] . 

In the direct numerical simulations of Eq. [21 r is incremented adiabatically in 
the course of the integration of the chain motion. Once [sup^, \xk\\ is inferior to 
a certain precision (i.e., 10~^) the external force is incremented. Before each 
increment, the chain configuration is recorded and once the chain has run over 
a distance L^., the integration is stopped. The latest anchored configuration 
corresponds to the strongest one and the associated external force is denoted as 
Tc, i.e., the static depinning threshold. We performed this type of simulations 
for different lattice aspect ratios, varying and Ly and for different obstacle 
densities ranging from 1 to 50 %. Various algorithms for the random numbers 
generator needed to build the obstacle array were tested and no significant 
difference was noticed in the end results. 
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Fig. 2. (color online) Strongest pinning configuration on a random array of obstacles (circles) for a chain 
length Ly = 2000 and a drag distance Lx = 500. In (a) the interaction cutoff is w = 0.5, the obstacle 
strength /jv/ = 0.005 and the obstacle density Cs = 7 %. In (b) w = 1, fj^j = 0.001 and Cs = 9 %. In the 
insets some segments are magnified with chain nodes marked as triangles (not visible on the main graphics). 
The X and Y axis scales differ for convenience of the plot. 



In Fig. [H we report the strongest chain configuration, obtained from the nu- 
merical simulations for a pinning strength Jm = 0.1. The critical profile is 
found to be wavy and to cross at least 40 lattice rows. In Figs. [2] (a) and (b), 
the critical chain profile is shown for smaller values of Jm, i-e., two orders 
of magnitude smaller than the one used in Fig. [1] We note that the entire 
string length is bounded by only two rows. The simulations evidence actually 
a well known feature for pinning of extended defects, namely weaker the ob- 
stacles fiatter the shape of the critical configuration. A perfectly rigid string 
would even experience a null force since then Jm would be negligible com- 
pared to the spring tension. However, as soon as some elasticity enters into 
play, the pinning strength becomes positive. The present work is essentially 
concerned with cases like those presented in Figs. [2] (a) and (b), where the 
elastic string shape is quasi-straight. In such situations, the string roughness 
is inferior or of the order of the inter-atomic spacing. The result shown in Fig. 
[T]only served us for comparison in order to introduce our problem. Such a case 
of wavy critical profile has been studied extensively, both through numerical 
simulations |32f33f34ll35ll36j and analytical works Pl6|7f8] . 



In the insets shown in Figs. [2] (a) and (b), it is worth noticing that along the 
rows that bound the spring chain, some holes appear in the obstacle distribu- 
tions. Hereafter, we dubbed such holes vacant site clusters. The sampling of 
such density fiuctuations along lattice rows plays a key role in the determina- 
tion of the critical drag force. 
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3 Vacant site cluster sampling theory 

3.1 The tightly bound chain 

In Fig. [2] (a), we noticed tliat for some parameters tlie strongest configuration 
of tlie spring chain remains tightly bound to the single lattice row at the back 
of the chain, i.e., most of the chain nodes are closer from the back row than 
from the next nearest one and thereby the chain does not cross several lattice 
rows. In the present section, a theory is devised to compute the critical drag 
force corresponding to critical configurations like the one shown in Fig. [2] (a). 
Further, the theory will be extended to cases with broader interaction cutoff. 

When the chain is tightly bound to the back row, the string can be viewed 
as quasi-straight, notwithstanding the bulges formed between rows. When 
w < 0.5, we can assume that the chain interacts with rows one by one and it 
is natural to work on the hypothesis that for such a system the strength of the 
random lattice is fixed by its more crowded row. To translate such a remark 
into some algebra, one needs to study the sampling of obstacles on a finite size 
lattice Lx x Ly. We notice that the purely random planar distribution follows 
Bernoulli's binomial law and the number of obstacles No involved into a single 
row of length Ly is then a random variable which probability is given by: 

pm = <ic^° (i-c.^-^°, (3) 

where C^^ = Ly\/No\{Ly — No)\. Such a statistical distribution can be approx- 
imated with a Poisson law in the limit of large Ly. However such a rounding 
yields some error for smallest Ly we are concerned with, so we keep the bino- 
mial formulation of Eq. [3l The probability for a row to involve less than 
obstacles is J2no<n Pi^o) and therefore in a set of L^ rows the probability for 
having a row with obstacles and (L^, — 1) rows with a number of obstacles 
inferior to is : 

(]{NJ = L, [p{Nj] [ Y: p{No)f^-'. (4) 

No<Nm 

The maximum number Nm fixing the number of obstacles in the denser row 
depends only on the lattice dimensions in each direction of space and on the 
overall obstacle density c^. It is easily computed numerically, paying attention 
to avoid overflows in factorials computation. The mean density in the denser 
row is then = N^/Ly. When an excess of vacant sites emerges at some 
place along the denser row, such segment is weaker than others where the 
obstacles are more crowded. Thence the spring chain should start the crossing 
at the largest vacant site clusters (VSC). The typical size of such VSC must 
now be determined. Actually the mean number of VSC in a row which obstacle 
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Fig. 3. (color online) Schematic representation of the model used for a quasi-straight spring chain tightly 
bound to a single lattice row. The small full circles represent the lattice sites, the large open circles represent 
the obstacles and the triangles arc for the spring chain nodes. The average number riv of vacant sites involved 
in the largest vacant site cluster is determined through Eq. [6] and Eq. [5] The average spacing m between 
the obstacles on both sides of the largest vacant site cluster is fixed by Eq. [7] 



The normalized probability 
- Cm)" while the probability 



density is fixed to Cm is Lysc = {cmLy — 1) ~ Nm 
to find a VSC with exactly n vacant sites is Cm(l 
for a VSC which size is inferior to n is [1 — (1 — Cm)"]- The probability to find 
a VSC of size n and {Lysc ~ 1) VSC with size inferior to n is proportional to: 



7(n) = LvsC [Cm(l - Cm)"] [1 - (1 - Cm,)"] 



(5) 



The mean size of the largest VSC in the denser row is thus: 



n. 



(6) 



Such a maximum VSC is surrounded by other VSC's that mean size is given 
by: [En<n, '^Cm(l-Cm)"]/[En<n, Cm(l-Cm)"] which for Convenience is denoted 
as (m — 1) with: 



1 (1 - Cm)"" 

m = riy z : — 

c 1 — fl — C 



(7) 



To compute the external force associated with the strongest binding row we 
consider the segment of vacant sites as embedded into a regular lattice 
of obstacles spaced by a mean distance m. Such a mean-field construction is 
illustrated within Fig. [3l where spring chain's nodes (triangles) are bound to 
the lattice sites occupied by the obstacles (large open circles). The array of 
obstacles is assumed to be centro-symmetric, so we ascribe the label to the 
center of symmetry which corresponds to the top of the bulge. We also define 
a new variable n = (1 + n^)/2 for convenience of the notations. Under the 
external applied force r, The force balance sheet, for say the left hand side of 
the chain leads to the set of equations: 
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Fo = -T - 2(xi - xo) 

Fi = -T - {x2 + a;o - 2xi) 

F2 = -T - (X3 + Xi - 2X2) 



n+m) 



Fn-l = -r - (X„ + Xn-2 " 2x„_i) 

Fn = -T - f{Xn) - {Xn+1 + Xn-l - 2x„) 
= -r - {Xn+2 + Xn- 2Xn+l) 

Pn+m—1 T~ {Xn+m ~\~ •^n+m— 2 2Xn+m—l) 

Fn+m fiXn+m) i^n+m+l ~l~ -^n+m— 1 2"^' 

-fri+m+1 (•^n+m+2 ~l~ X^+m 2x^_|_m+l)) (8) 

and in principle the series of equations repeats up to the chain boundaries 
with increment of subscripts. We assume that the mechanical equilibrium is 
satisfied for all nodes j situated in between obstacles. Then we have Fj = 
but for j G [n, n + m, ..n + pm]. For the segment j G [0, n], it is easy to show 
by recurrence that : xj — xq = —Tp/2. Thence the chain shape is parabolic 
between n and —n. For j G [n, n+m], we proceed the same and find x„ = 

—rj{n + j/2) — [Fn + f{xn)]j which fixes the segment end to 

Xn+m = Xn- Tm{n + m/2) - m[Fn + f{Xn)]- (9) 
The same can be iterated once again for j G [n + m, n + 2m] which leads to 

Xn+2m - Xn+m = -Tm{n + 3m/2) - (F„ + Fn+m + f{Xn) + f{Xn+m))m. The 

set of equation on the positions Xn+jm is then: 



Xn+m = Xn - Tm{n + m/2) - f\xn)m 
Xn+2m = Xn+m ~ Tm{n + 3m/2) - [F„ + 

Pn+m + f{Xn) + f{Xn+m)]m 



Xn+pm — Xn+{p-l)m ~ Tm{n + {2p — 1)"^) 

Subtracting the two latest equations yields: 

P ^mXn+pm n/ \ (11\ 

^ n+pm Tm J [Xn+pm ) \^^) 

m 

where ApXn+pm = (a;„,+(p+i)m + Xn+{p-.i)m - "^Xn+pm) IS the discrete Laplacian 
applied to the p subscript. When the entire chain is at mechanical equilibrium 
Fn+pm = for all p. Far enough from the n^-VSC (i.e., the VSC with riy vacant 
sites), the solution for Xn+pm tends asymptotically to a constant Xoo such as 



- m^[Fn+jm + f{Xn+jm)]i^^) 
7=0 
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rm = — /(xoo) and therefore: 

2 arccos (-rm//M) , 47r 
3^00 = ^ cos( + —). (12) 

We can expend linearly Eq.[TT]for the far enough sites such as the displacement 
Xn+pm writes Xn+pm = x^ + Cp and f{xn+pm) = f{xoo) + f'{xoo)ep. Then, at the 
equihbrium Eq. [TT]yields [Apep = —f'{xoo)'mep] and thence ep is an exponential 
function: [ep = eoexp(— ap)] which the exponent a verifies 



a = ±2 ash(y'-/'(xoo)m/2). (13) 

Since the chain displacement is bounded, we are solely concerned with so- 
lutions such as (ap) > 0. The sum of the whole set of equations in Eq. [TOl 
provides another relation between r and the nodes position 

Xn+pm.} the 

condition that Fn+pm = for all p: 

p-i 

Xn - Xn+pm = m[Tp{n + pm/2) + "^{p ~ j)f{Xn+jm)], (14) 

i=o 

which after expanding f{xn+jm) as a Taylor series around Xoo and keeping 
only the terms linear in p provides us with an equation which relates r to eo: 

"^o r f'{Xoo) , ^of"{xoo) I 4f"'{xoo) , , 



(n-m/2)'(l-e-") 2(1 - e'^") 6(1-6" 



The critical chain configuration is reached when the Hessian associated with 
Eq. [11] has a singular eigenvalue. This allows us to determine the critical 
value for eo. Actually we found that finding the Hessian singular eigenvalue 
is equivalent to find the maximum of Eq. [15] for r with respect to eo. The 
solution for the critical bulge is then: 

^ / (3^00) \j f '{Xoo)"^ 4/"'(Xoo)/'(3^oo) 3(]^„e-3a^(]^lg-~ 

'° = "^^^ " '"'"^ 2r(xoo)(l-e-2") • 

(16) 

Combining the solutions for Eqs. [12], [131 IE] and [16] allows us to determine 
the maximum pinning force associated with N^, the number of obstacle in 
the denser row. For this reason, we denote such a maximum as r(iVm). The 
set of equations giving r(A^m) can be solved recursively. Starting with a small 
enough trial solution for r = ro, we compute the corresponding quantities Xoo 
and a from Eq. [T2] and from Eq. [T3] Then eg is derived from Eq. [TH] and the 
corresponding value of r from Eq. [15] If the so obtained quantity is larger than 
the initial value Tq then the latter is incremented and we proceed the same up 
to find identical values for r and To. The end result gives the required T{Nm) 
to a precision fixed by the trial solution increment. 
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Fig. 4. (color online) External force required to drag a spring chain of length Ly over a distance Lx, for 
a pinning strength /j;/ = 0.005 and a cutoff w = 0.1 in (a) and w = 0.5 in (b). The symbols represent 
the data obtained through the simulations described in Sec. [21 for different lattices (see figures legend). 
The continuous lines correspond to the predictions made through the theory detailed in Sec. \3\ for same 
parameters as those used in simulations. Colors of symbols and lines correspond one to one. 
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Fig. 5. (color online) External force required to drag a spring chain of length Ly = 1000 along a distance 
Lx = 100, for different obstacle pinning strengths (see legend) and the same interaction range w = 0.5. 
The different symbols represent the simulations data and the continuous lines correspond to the prediction 
made through the analytical theory detailed in Sec. |3] for same parameters. 

The critical pinning force Tc of the random lattice is approximated by averaging 
T{Nm) over A^^: 

n = Y.P{NmMN^), (17) 

where (3{Nm) has been given in Eq. HI The previous theory is compared to sim- 
ulations data in Figs. H] (a) and (b) and in Fig.[5]for different lattice dimensions, 
different pinning forces, varying w and J'm- A quantitative agreement has been 
obtained between theory and simulations, although no adjustable parameters 
are involved. According to same type of comparisons but for larger /m values, 
the previous analytical work proves relevant for Jm smaller than roughly 0.03. 
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position of the chain / lattice X-axis 

Fig. 6. (color online) Strongest pinning configuration for fj^j = 0.03, Ly = fOOO, = fOO, a density 
Cs = 9 % and an interaction cutoff w = 0.5. The inset shows a chain segment the nodes of which are marlced 
with triangles. Only obstacles close from the string have been reported as circles. 

The theory predictions worsen for cases where the critical configuration crosses 
few lattice rows. As an example, for = 0.03 in Fig. [5l the VSC theory is 
found to become less accurate for low densities. In such conditions, the string 
profile at the depinning transition corresponds to a kinked shape, as shown 
in Fig. [6l different from the wavy profile shown in Fig. [1] and from the quasi- 
straight ones shown in Figs. [2] (a) and (b). From Fig. it can be seen that the 
discrepancy increases as the density decreases while the theoretical predictions 
remain accurate for more concentrated obstacle distributions. As /m increases 
above 0.03, the deviation between theoretical predictions and simulations data 
is shifted toward higher densities. In this range of parameters, the system 
undergoes a bifurcation, not treated in the present work. 

Other authors [32][33] noticed that the maximum pinning force of a random lat- 
tice was dependant of the drag distance. Concerning the weak pinning points 
studied here, it is thus of some interest to explore the variation of the critical 
drag force with lattice dimensions. The critical drag force was found to vary 
proportionally to [ln(L2;)]""= where varies with all parameters. For instance 
when w = 0.5, /a/ = 0.01 and Ly = 1000, we found = [0.45 - 0.036 ln(c,)]. 
One can thus conclude that the dependence is very weak since a frac- 
tional power of a logarithm is a rather wise function. The maximum drag 
force depends not only on but also on Ly. The critical drag force varies 
proportionally to a constant plus the function [ln(Ly)]~"!' where ay decreases 
when Cs increases. For instance when w = 0.5, /m = 0.01 and L^,. = 4000, we 
found ay = [0.26 — 0.631n(cs)]. According to direct simulations, the Tc depen- 
dence in Ly seems to weaken when J'm becomes large enough for the critical 
configuration to differ from the quasi-straight line as in Fig. [6] or Fig. [H 

For fixed lattice dimensions, the adjustment of a density power law for the 
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theoretical critical drag force allows to establish some comparison with the 
SSH analytical theory. The density power law fit is found very close from a 
linear variation and may even overpass slightly the unitary exponent in some 
situations depending on the lattice geometry. For instance, from Fig. |l](b), we 
worked out by a curve fitting with the form Tc oc (c^), an exponent t] = 1.13 
for Ly = 4000 and = 25 whereas for Ly = 500 and = 1000, t] = 0.946 
was obtained. The effective density exponent rj is therefore size dependent. 
A density exponent close from unity corresponds to the FMS theory for the 
hardening of concentrated solid solutions. The CRSS linear dependence in the 
solute concentration was also noted in atomistic calculations on the model 
solid solution Ni(Al) [28]. 

3.2 Extension to broader cutoff' 

To extend the theory developed previously for a short cutoff w, we first con- 
sider the case w ~ 1, which allows us to treat the spring chain interaction 
with only two rows: the row at the back and the next nearest one. Such a 
critical configuration can be seen from direct simulations as reported in Fig. 
H] (b). To compute Tc, we still employ the same model as depicted in Fig. [31 
However, the interaction with the next nearest row may weaken the pinning 
on the back row. One ascribes to each site j of the next nearest row, a force 
g{xj) = Cvf{xj — 1) where c„ is the density of obstacles in the next nearest row. 
For the situation shown in Fig. [31 the force balance sheet writes as follows: 

Fo = -T - g{xo) - 2(xi - a;o) 

Fi = -T- g{xi) - {x2 + Xo- 2xi) 

F2 = -T - g{x2) - (xs + xi- 2x2) 

Fn-l = -T - g{Xn~l) - {Xn + Xn-2 " 2x„_i) 

Fn = -T - g{Xn) - f{Xn) " {Xn+1 + X^-l - 2x„) 
Fn+1 = -T - g{Xn+l) - {Xn+2 + X^ - 2x„+i) 

Fn-i-m—1 fi'('^n+m— 1) (-^n+ni ~l~ ■^n+m— 2 2Xn+m—l) 

F-n+m T~ g{Xn+m) f^Xn+m) (-^n+m+l ~l~ Xn+m~l 2Xn+m) 

Fn+m+1 = ~T ~ g{Xn+m+l) ~ {Xn+m+2 + Xn+m ~ 2Xn+m+l) 

(18) 

To estimate the force exerted by the next nearest row upon the spring chain, 
the function g{xj) is approximated with a step function: g{xj) = g{xo) if 
j E [0;n - (m - l)/2[ and g{xj) = g{xn+pm) if j e [n + pm - {m - l)/2;n + 
pm + (m+ l)/2[. We then assume that the mechanical equilibrium is achieved 
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for all nodes excepted those aligned with some obstacles of the back row. This 
leads then to the following set of equations: 



T = -g{xo) - 2(xi - Xq) 

T = -g{xo) - (x2 + a;o - 2xi) 

r = -g{Xn) - {Xn + Xn-2 - 2x„_i) 

T f{-^n) (-^n+l ~^ -^n— 1 Sx^) 

r = -g{Xn) - {Xn+2 +Xn- 2x„+i) 



n+m— 2 



- 2x 



n+m—l 



-^n+m fiS^n+m) diS^n+m) ~l~ •^n+m—l '^•^n+rri) 

T = —g{Xn+m) ~ {Xn+m+2 + Xn+m ~ 2x„+m+l), (19) 

and the equations repeat up to the chain boundaries by incrementing sub- 
scripts. By applying recurrence, it is possible to reduce the previous set of 
equations to a smaller one, concerning only the regular array of obstacles in 
the back row: 



p-i 

Xn+pm = Xn+{p-i)m " rm{n + {2p - l)m/2) - m\Y^\Fn+jm + 



p-2 

+m'^g{xn+jm) + g{xo){n - m/2) 

j=0 

-m giXn+{p-l)m)- 



— 1 



g (^^n+pm) 



(20) 



Subtracting the equation for rank (p — 1) from the one at rank (p) gives: 



n+pm ■ 



^m-^n+pm TTfl n^'fi^^n+pm) 



m 



9{Xn+{p+l)m) + 9{Xn+{p-l)m) 



3m2 + 1 



9i.Xn+pm), (21) 



while for p = 0: 



Fn = Xn - Xn+m " Tm{n + m/2) - m[f{xn) + 9{xo){n - m/2)] 



rn? — 1 



8 



-g[Xn+ni) - m- 



3m — 1 



-9[Xn) 



(22) 



When the equilibrium is achieved Fn+pm = for all p, which allows us to 
;he asymptotical solution of E 

verifies: 



deduce the asymptotical solution of Eq. [2T] as a constant Xn+pm which 



rm = -/(xoo) - mg{xoo)- 



(23) 
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Fig. 7. (color online) Same as in Fig|3]but for a cutoff w = 1 and a pinning strength = 0.001. 

The solution to Eq. [23] is the positive real root of a third order polynomial 
equation: 



1 + mcv 
with Aa 

and Br^, 



cos arccos 



-5. 



27 X , 4:^T^ 



— 1 — mCt, — 



1 



+ 



1 + mCy 
3W^" 



^2(1 + mc^ 



w(l + mCy) 



(24) 



If (1 — Xoo) > w then we can set c„ = in the previous equation which 
leads to Eq. [121 valid for small w. We now expend linearly Eq. [21] around Xoo 
to determine an approximation of nodes position as Xn+pm = Xoo + e^. The 
equation on is: 



-m/'(xoo)ep 



3m^ 



— 1 



4 



g'{xoo) ep+i + ep_i . (25) 



The solution is an exponential function ep = eoexp(— ap) with a dispersion 
relation: 



/ 2 - mfixoo) - (3m^ + l)g'{xoo)/A ^ 
[ 2 + ^.9'(Xoo) ^' 



a = ±ach( - — Y^v'^-^^- ). (26) 

At this stage, it is of some interest to work out the maximum of the spring 
chain position. For the segment situated along the larger hole (see Fig. [3]), 
j G [0,n — (m — l)/2[, the set of very first equations in Eqs. [19] leads to: 



rn 



rrr — 1 



Xfi Xq 



-9[Xn) 



\n — m' 



74 + 1/4]. 



(27) 
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Then xq can be expressed as a function of x„ since xq is actually the positive 
root of a third order polynomial: 



Xo = l + 2w]^ 
with Aq 



I Aq + w 



cos 



arccos 



Co 



27AI 



2Ao\ iAo + 



w 



/3 + 47r/3 



n 



mV4 + 1/4 



and Co = Xn — I 



Tfi m 
^ + — 



'-9{Xn)- 



(28) 



We also need to express the first derivative of xo against Xn which according 
to Eq. EH gives: 



dxo 



1 + r^g'jxn) 

[^2 _ (^2 _ l)/4]c,'(a;o)/2' 



(29) 



The sum of the equations in Eq. [20] from rank 1 to rank p, taken sufficiently 
large, leads to an equation which the linear term in p is: 



r 



-9{xo) + 



[n — m/2] 



5m 



eo- 



(/'(a^oo) +mg'{x^)) 



(/"(Xoo) + mg'\x^)) 3 (/"'(Xoo) + mg"\x^)) 



2(1 



-2a^ 



6(1 



(30) 



The latter equation is similar to Eq. [15] obtained for the tightly bound chain, 
but includes the interaction with the next nearest row. The maximum r in 
Eq. [30] against eo corresponds to the critical strength which provides a tran- 
scendental equation on eo: 



[n-m/2\g'{xo) 



dxo 



if'ixoo) + mg'{xoo)) (/"(a^oo) + rng'^Xoo)) 
eo — 



1 - e-" 
,{f"{x^)+mg"'{x^)) 



l-e 



-e, 







2(1 



-3a\ 



(31) 



where Xn = x^o + Co; and dxo/dx^ are given in Eq. [23 and Eq. [23 This 
complete our computation of the maximum pinning force associated with 
and Cy. The corresponding value of r is therefore related to A''^ through Eq. [6] 
as well as to the number of obstacles iV„ = CvLy. The maximum force is now 
considered as a function of both Nm and A''^ and it is denoted by T{Nm, N^)- 
The probability to find a couple of rows which actually consists of a back row 
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Fig. 8. (color online) External force required to drag a spring chain of length Ly = 200 over a distance 
Lx = 1500, for different obstacle pinning strengths (see legend) and the same interaction range w = 1. 
The different symbols represent the simulations data and the continuous lines correspond to the predictions 
established through the analytical theory detailed in the text. 

with Nm obstacles and a front row with A'^ obstacles as being the strongest 
configuration among rows is written as: 

^(iV^,iV,) = [p{Njp{N,)] [1 - E p(iVi)p(iV2)]^^-\ (32) 

r{NuN2)>TiN^,N.,) 

where the function p{N) is the binomial given in Eq. [31 Thence the average 
critical depinning is given by: 

Tc= r{N„^,N,)e{N„,,N,). (33) 

The same statistical treatment can be applied to the case w < 0.5 and yields 
the same results as presented in Sec. 13.11 since the critical force given in Eq. 
[15] is then independent of A''^,. 

In Fig. [7] and Fig. [HI the theoretical predictions are compared to simulations 
data. It is worth noticing that the critical force variation against density cor- 
responds to a more convex curve than for shorter cutoff w < 0.5. An analysis 
in term of power law fit leads to an exponent smaller than unity. For instance 
in Fig. [7], for Ly = 500 and = 1000, the effective density exponent is around 
a = 0.65. The effective density exponent of Tc is found to decrease when w 
increases (see further Sec. 13. 3p . In Fig. [7], the computations have been per- 
formed for different aspect ratios of the lattice and show the same trend as 
in Fig. [H for a shorter w, the critical force density exponent increases with 
Lj. and decreases with Ly. Although a satisfactory agreement is obtained for 
the different lattices, we notice that the pinning force Jm is smaller than in 
Fig. [H Actually the field of validity for the VSC theory is narrower for w = 1 
than for w < 0.5. The limit of application for the theory decreases as the 
interaction cutoff increases. Such a limit also varies with the lattice dimen- 
sions: it decreases when Ly increases and when L^ decreases. For instance. 
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with Ly = 200 and = 1500, the theory proves to be efficient as seen from 
Fig. [8] where the comparison has been performed for different pinning forces 
the maximum of which is J'm = 0.006. The same computations carried out for 
Ly = 2000 and = 500, keeping constant both w and /m yield much worse 
results, in particular for the low densities. This can be understood comparing 
the critical proffies for both geometries. Whereas the chain profile is quasi- 
straight for the former the critical profile is wavy for the latter, i.e., similar to 
the one seen in Fig. [TJ As for short cutoffs, the system undergoes a bifurcation 
passing from a quasi-straight critical proffie to a wandering one. The change 
in critical proffie occurs for smaller J'm with w = 1 than for w < 0.5. 



3.3 Extension to intermediary and still broader cutoffs 



In previous computations, we developed the pinning force function as a Taylor 
series around Xoo, i.e., the asymptotic solution for the chain nodes position. 
For intermediary cutoff, namely 0.5 < w < 1, we must be aware that such a 
development cannot be used to approach the non-analytic force function since 
the cutoff occurs right in between two rows. Then Xoo is assumed to remain far 
from the cutoff interaction with the next nearest row. The asymptotic position 
Xoo is thus given by Eq. [12] and a by Eq. [131 The top of the bulge is assumed 
to be situated above the cutoff distance from the next nearest row. Then the 
equation on xq is same as in Sec. 13.21 and can be solved analytically as shown 
in Eq. [28l The equation relating r to cq must be rederived. The sum of the 
equations in Eq. [20] from rank 1 to rank p, taken sufficiently large, leads to an 
equation which the linear order in p is now: 



_ _ / \ ^ \ f (^oo) 2 / (^oo) , 3 / i^oo) 

+ ^9{Xn+jm)\: (34) 

i<ib 

where is such as {xn+it,m > I — w). The quantity lb corresponds to the 
segment of the bulge which overpasses the interaction cutoff with next nearest 
row. Such equation for r holds provided that lb remains small compared to 
Ly/m. 

In Fig. [9] the predictions from the present development are compared to sim- 
ulations data for w = 0.7. We reported the results obtained from simulations 
with different obstacle densities. For some densities, we carried out the simula- 
tions with different random distributions. The simulations data are distributed 
around the theoretical predictions which con&ms the robustness of the theory. 

The previous computations can also be extended to an interaction cutoff larger 
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Fig. 9. (color online) External force required to drag a spring chain of length Ly = 1000 over a distance 
Lx = 200, for a pinning strength /m = 0.001 and different interaction range ui (see legend). The different 
symbols represent the simulations data and the continuous lines correspond to the prediction made through 
the theory detailed in the text. 

than the unit cell parameter (i.e., w > 1). In the theory, the function g{xj) = 
Cvf{xj — 1) is substituted with g{xj) = c^f{xj — 1) + h{xj) where h{xj) = 
Csf{xj — 2) accounts for the pinning force due to the still further row, i.e., 
the second next nearest one. The top of the bulge may interact with such a 
row meanwhile the back of the chain is assumed to remain far from it. Then 
the equations found in Sec. 13.21 for Xoo and a still hold and xq is again given 
by Eq. [27| but substituting the function g{x) with the expression g{xj) = 
Cvf{xj — 1) + h{xj). The equation on r follows: 




r 



-g{xo) + 



\n 



(m + 11 



-g[Xc 



eo- 



— e 



5m 1 — e 

2 U"i.Xoo) +mg"{xoo)) 
-0 2(1 - e-2«) 
,{f"'{x^)+mg"'{x^)) 







6(1 -e- 



■3a-] 



EM 



n+jm ) 



(35) 



3<h 



where corresponds to the top of the bulge overpassing the distance cutoff 
with the second next nearest row, defined as [xn+i^m > 2 — w). In Fig. [9l 
the theoretical computation for against the obstacle density is compared 
with simulation data for {w = 1.1) and for different random distributions. 
Although the theory could certainly be improved by accounting for the effect 
of the nearest row at the back of the chain, the agreement with simulations 
data proves quantitative. 

The results obtained in both cases, w = 0.7 and w = 1.1 show that Tc decreases 
with w. This was confirmed by other simulations performed for still larger w, 
up to w = 2. Such a decreases is opposite to the one predicted within the 
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MNL theory for SSH. Moreover according to the MNL theory, the density 
exponent of Tc is fixed to t] = 2/3, whereas according to the VSC samphng 
theory such an exponent depends on lattice dimensions (see Sec. 13.21) and w, 
as welL For instance, in Fig. [HI a power law fit on the theoretical predictions 
yields and effective exponent 77 = 0.945 for w = 0.7 while for u; = 1.1 we 
obtained r] = 0.827. The variations 77 with lattice dimensions, along with 
its value larger than 2/3 and its w dependence demonstrate that the MNL 
theory does not apply to the systems studied here. The same conclusion holds 
for other SSH theories predicting a constant density exponent of the critical 
depinning threshold. 



4 Summary and perspectives 

In the present study, we addressed a discrete version of a paradigmatic prob- 
lem, namely the depinning threshold of an elastic string on a random substrate. 
For a planar distribution of weak pinning points, a theory was devised to com- 
pute the applied force required to drag the one-dimensional elastic manyfold 
over a disordered potential landscape with various aspect ratios. The the- 
oretical predictions were found accurate provided the critical configuration 
remains close from a quasi-straight line bounded between two lattice rows. 
The strongest pinning configuration was shown to hinge on the denser lattice 
rows in which the largest vacant site clusters (VSC) are bounded in size. Such 
maximum VSC correspond to the weakest lattice defects on which the critical 
depinning proceeds. The mean size of the critical VSC is determined through 
an expression involving only lattice dimensions and the overall planar obstacle 
density Cg. The theory allowed us to account for the finite lattice size effects 
in a quantitative manner. The typical variations of the critical applied force 
against the chain length and the drag distance both yield a logarithmic power 
law. For a fixed lattice geometry, and an interaction cutoff inferior to half the 
lattice cell parameter the pinning strength was found close to being linear in 
Cg. The effective density exponent rj of the depinning threshold was found to 
depend on the lattice geometry and to decrease as the interaction cutoff w 
increases. 

In some atomistic studies bearing on dislocation in a model Ni(Al) solid solu- 
tion, i.e., with a rather marked size effect [2H], the maximum force exerted by 
isolated solute atoms on a dislocation segment was found of the order of 0.04 
nano-Newton (nN). The elastic contribution to the line tension of a screw dis- 
location, computed within the isotropic elastic theory [37] is F = rsln{R/b) 
where b is the Burgers vector, R is the inter-dislocation spacing and where 
the pre- logarithmic factor is given by F^ = /i6^(l + z/)/47r(l — z/), with /i as 
the shear modulus and u as the Poisson coefficient. The distance R is related 
to the dislocation density pd hj R = l/^/pd which for a standard density in 
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deformed metals pd = 10^^ m~^, gives R = 1 fim. Then, with fiNi = 74600 
MPa and u = 0.28, the screw dislocation line tension in Ni is F = 6.1 nN 
which gives a ratio between the obstacle pinning strength and the line tension 
of Im = 0.007, of the order of /m studied in the present work. A solute atom 
dilation smaller than the one for Al in a Ni matrix or a smaller dislocation 
density could even yield smaller Jm- In regard of the simplicity of the elastic 
manyfold model, the present study requires to be extended to a more realistic 
model involving some of the atomistic details as for instance: (i) a position 
dependant obstacle strength, (ii) a dissociated core, and (iii) mixing of differ- 
ent types of obstacles. Such extensions would allow to address longstanding 
issues in SSH theory. 
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Depinning of a discrete elastic string from a 
two dimensional random array of weak 
pinning points 
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Abstract 

The present work is essentially concerned with the development of statistical theory 
for the low temperature dislocation glide in concentrated solid solutions where atom- 
sized obstacles impede plastic flow. In connection with such a problem, we compute 
analytically the external force required to drag an elastic string along a discrete two- 
dimensional square lattice, where some obstacles have been randomly distributed. 
The corresponding numerical simulations allow us to demonstrate a remarkable 
agreement between simulations and theory for an obstacle density ranging from 1 
to 50 % and for lattices with different aspect ratios. The theory proves efficient on 
the condition that the obstacle-chain interaction remains sufficiently weak compared 
to the string stiffness. 

Key words: depinning transition, dislocation, solid solution hardening 
PACS: 61.72.Lk,74.25.Qt,64.60.An 



1 Prom the solid solution strengthening theory 



The statistical theory for solid solution hardening (SSH) emerged from the 
seminal works of Sir N. Mott[?] and his near colleagues, F.R.N. Nabarro[?,?] 
and J. Friedel [?]. The early analytical theory, perfected and extended by other 
contributors, as for instance R. Fleischer, R. Labusch and T. Suzuki [?,?,?,?,?] 
applies to substitutional alloys where the solute atoms can be considered as 
immobile during the dislocation glide, by contrast to the cases where dislo- 
cations may drag along an atmosphere of fast diffusing impurities. In face 
centered cubic (fee) alloys, the critical resolved shear stress (CRSS) was then 
expected to increase in proportion to with Cs as the atomic concentration 
of solute atoms and r] as an exponent depending on the assumptions made on 
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the interaction between dislocations and foreign atoms: 77 = 1/2 in Friedel- 
Fleischer (FF) theory while 77 = 2/3 in Mott-Nabarro-Labusch (MNL) theory 
and ?7 = 1 in Friedel-Mott- Suzuki (FMS) [?,?]• Within analytical theory for 
SSH, the dislocation is thought of as a continuous elastic string impinged on 
a two-dimensional (2D) random static potential. The depinning transition in 
such a model is a typical issue of statistical physics, belonging to a broad 
class of problems concerned with extended interfaces motion in heterogeneous 
materials!?,?,?,?,?,?,?,?,?,?,?,?,?]. 

The recent developments of three-dimensional atomistic simulations (3D- AS) 
allowed to work on more realistic models for dislocations in solid solutions [?,?,?, ?,?,?,?,?]. 
Though 3D-AS confirmed that a large part of the dislocation pinning hinges 
on the impurities situated in the crystal planes that bounds the dislocation 
glide plane[?], the simulations revealed also the complexity of the dislocation- 
obstacle interaction. In fee alloys, the geometry of the dislocation core, dis- 
sociated in two Shockley partials separated by a (111) stacking fault ribbon 
undermines the simple picture of an elastic line in interaction with a single 
type of obstacles, as stated in the basic version of SSH theory. Instead, the 
pinning forces differ according to partials and to the obstacle positions, i.e., 
above or below the glide plane [?,?]. 

On the other hand, the nanometric scale of the atomistic simulations, a strin- 
gent limit imposed by the computational load, hinders the direct extrapola- 
tion of simulation results to macroscopic samples. A multi-scale approach is 
thence required to link the atomistic studies to the realm of materials sciences. 
A manner to proceed consists in incorporating some of the atomic ingredients 
revealed in 3D-AS to a discrete version of the elastic string model that remains 
tractable even for large dimensions [?]. The discretization of the string model 
should allow to transfer of the atomic details. 

In the present paper, the impact of the discretization on the depinning transi- 
tion is analyzed thoroughly. The elastic string is replaced with a discrete spring 
chain, the nodes of which move on a 2D square lattice and interact with some 
pinning points randomly distributed on lattice sites. This very simple model 
allows us to devise an analytical theory which accounts for the discreteness of 
the obstacle distribution and thus opens a promising perspective to integrate 
more of the atomic details. In order to demonstrate the accuracy of the theory, 
we compute directly the critical external force within numerical simulations 
applied to the discrete string model. Theory and simulations agree remark- 
ably well for a broad range of model parameters, e.g., (i) the in-plane obstacle 
density c^, (ii) the lattice size in every direction of space, (iii) the maximum 
pinning force /m and (iv) the potential interaction cutoff w that characterizes 
the obstacles. The theoretical predictions are proved reliable on the condition 
that /m and w remain smaller than certain bounds varying with c^. For a dense 
distribution of weak pinning points, the critical configuration of the chain is 
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found to be a quasi-straight line parallel to the atomic rows. The depinning is 
then shown to occur at some vacant site clusters (VSC) which the typical size 
is explicitly related to the lattice dimensions in both directions of space. The 
external force required to drag along the spring chain over a finite distance 
reflects such a size dependence. Noteworthily the effective density exponent 
rj is also found to vary with lattice dimensions, in contrast with expectations 
drawn on standard SSH theory. 

Our report is organized as follows. In Sec. II, the spring chain model is intro- 
duced and the direct numerical computations are described. In Sec. Ill, the 
statistical theory is derived and compared to the numerical data. The results 
are resumed and commented in Sec. IV. 



2 The phenomenological spring chain model 

The model proposed hereafter belongs to the wide class of elastic interface 
models, extensively studied in statistical physics [?,?,?,?,?,?,?]. A one dimen- 
sional elastic string is discretized with a spatial step b, equivalent to the short- 
est interatomic distance in solids. Each node of the discrete chain is bound to 
its first neighbor by an harmonic spring of strength F. The two quantities, b 
and r are chosen to scale distances and forces, respectively. The spring chain 
nodes move along the column of a square lattice. The size of the lattice in the 
direction of the chain is denoted as Ly whereas the distance over which the 
chain is dragged is L^. The 2D random array of obstacles is constructed by se- 
lecting the occupied lattice sites, up to a number of obstacle equals to CsL^Ly, 
where Cg is the obstacle density. Since the depinning process occurs when the 
chain nodes pass the force maximum, the interaction potential is expended as 
a polynomial function in the vicinity of such a maximum. Assuming that the 
interaction is attractive and that the potential is symmetric with respect to 
its minimum, we obtain a polynomial function of at least fourth order: 

V{x) = Vo{x^/w^ - 1)^ for |a;| < w 

V{x) =0 for \x\ > w, (1) 

which corresponds to a force f{x) = — 4Vo(x^/tf^ — l)x/w'^, with a maximum 
value /m = 8| Vo|/(3v^w), attained when x = ±w/y/3. The chain nodes inter- 
act solely with obstacles situated in the column along which they may glide. 
The polynomial force with a distance cutoff w is obviously very far from the 
dislocation-solute interaction, characterized by a decrease of Coulomb type. 
Hereby we describe only the local potential variation yielded when a solute 
atom visits a dislocation core. The parameter w fixes how the interaction de- 
creases in the vicinity of the force maximum. Both Jm and w can be extracted 
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position of the chain / lattice X-axis 



Fig. 1. (color online) Strongest pinning configuration of the spring chain on a random array of obstacles 
(circles) for fm = 0.1, Ly = 2000, Lx = 500, a density Cs = 16 % and an interaction cutoff w = 1. 
Only obstacles close from the chain have been reported for clarity. Inset shows a magnification of obstacles 
(circles) and nodes (triangles) of a chain segment. X and Y axis have different scaling for convenience of 
the plot. 

from atomistic data as those reported in [?,?]. 

The dimensionless over-damped Langevin dynamics for the chain node k is 
given by: 

X, = [xfc+i + x,.^ - 2x,] + r - E 4^^ (^i^iMl! _ 1), (2) 

where Xk is the position of the node k, r is the external force and Sk,i is the 
coordinate of the ith obstacle in the kth row. For the weak pinning forces 
we are concerned with, the chain strain remains very small such that the 
anharmonic terms in the spring tension have been neglected. Properly scaled, 
the continuous version of the spring chain model served in the development of 
the SSH theory [?,?,?,?]. 

In the direct numerical simulations of Eq. [21 r is incremented adiabatically in 
the course of the integration of the chain motion. Once [sup^, |xfc|] is inferior to 
a certain precision (i.e., 10~^) the external force is incremented. Before each 
increment, the chain configuration is recorded and once the chain has run over 
a distance L^, the integration is stopped. The latest anchored configuration 
corresponds to the strongest one and the associated external force is denoted as 
Tc, i.e., the static depinning threshold. We performed this type of simulations 
for different lattice aspect ratios, varying and Ly and for different obstacle 
densities ranging from 1 to 50 %. Various algorithms for the random numbers 
generator needed to build the obstacle array were tested and no significant 
difference was noticed in the end results. 
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Fig. 2. (color online) Strongest pinning configuration on a random array of obstacles (circles) for a chain 
length Ly = 2000 and a drag distance Lx = 500. In (a) the interaction cutoff is w = 0.5, the obstacle 
strength /jv/ = 0.005 and the obstacle density Cs = 7 %. In (b) w = 1, fj^j = 0.001 and Cs = 9 %. In the 
insets some segments are magnified with chain nodes marked as triangles (not visible on the main graphics). 
The X and Y axis scales differ for convenience of the plot. 



In Fig. [H we report the strongest chain configuration, obtained from the nu- 
merical simulations for a pinning strength Jm = 0.1. The critical profile is 
found to be wavy and to cross at least 40 lattice rows. In Figs. [2] (a) and (b), 
the critical chain profile is shown for smaller values of Jm, i-e., two orders 
of magnitude smaller than the one used in Fig. [1] We note that the entire 
string length is bounded by only two rows. The simulations evidence actually 
a well known feature for pinning of extended defects, namely weaker the ob- 
stacles fiatter the shape of the critical configuration. A perfectly rigid string 
would even experience a null force since then Jm would be negligible com- 
pared to the spring tension. However, as soon as some elasticity enters into 
play, the pinning strength becomes positive. The present work is essentially 
concerned with cases like those presented in Figs. [2] (a) and (b), where the 
elastic string shape is quasi-straight. In such situations, the string roughness 
is inferior or of the order of the inter-atomic spacing. The result shown in Fig. 
[T]only served us for comparison in order to introduce our problem. Such a case 
of wavy critical profile has been studied extensively, both through numerical 
simulations [?,?,?,?,?] and analytical works [?,?,?,?]. 



In the insets shown in Figs. [2] (a) and (b), it is worth noticing that along the 
rows that bound the spring chain, some holes appear in the obstacle distribu- 
tions. Hereafter, we dubbed such holes vacant site clusters. The sampling of 
such density fiuctuations along lattice rows plays a key role in the determina- 
tion of the critical drag force. 
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3 Vacant site cluster sampling theory 

3.1 The tightly bound chain 

In Fig. [2] (a), we noticed tliat for some parameters tlie strongest configuration 
of tlie spring chain remains tightly bound to the single lattice row at the back 
of the chain, i.e., most of the chain nodes are closer from the back row than 
from the next nearest one and thereby the chain does not cross several lattice 
rows. In the present section, a theory is devised to compute the critical drag 
force corresponding to critical configurations like the one shown in Fig. [2] (a). 
Further, the theory will be extended to cases with broader interaction cutoff. 

When the chain is tightly bound to the back row, the string can be viewed 
as quasi-straight, notwithstanding the bulges formed between rows. When 
w < 0.5, we can assume that the chain interacts with rows one by one and it 
is natural to work on the hypothesis that for such a system the strength of the 
random lattice is fixed by its more crowded row. To translate such a remark 
into some algebra, one needs to study the sampling of obstacles on a finite size 
lattice Lx x Ly. We notice that the purely random planar distribution follows 
Bernoulli's binomial law and the number of obstacles No involved into a single 
row of length Ly is then a random variable which probability is given by: 

pm = <ic^° (i-c.^-^°, (3) 

where C^^ = Ly\/No\{Ly — No)\. Such a statistical distribution can be approx- 
imated with a Poisson law in the limit of large Ly. However such a rounding 
yields some error for smallest Ly we are concerned with, so we keep the bino- 
mial formulation of Eq. [3l The probability for a row to involve less than 
obstacles is J2no<n Pi^o) and therefore in a set of L^ rows the probability for 
having a row with obstacles and (L^, — 1) rows with a number of obstacles 
inferior to is : 

(]{NJ = L, [p{Nj] [ Y: p{No)f^-'. (4) 

No<Nm 

The maximum number Nm fixing the number of obstacles in the denser row 
depends only on the lattice dimensions in each direction of space and on the 
overall obstacle density c^. It is easily computed numerically, paying attention 
to avoid overflows in factorials computation. The mean density in the denser 
row is then = N^/Ly. When an excess of vacant sites emerges at some 
place along the denser row, such segment is weaker than others where the 
obstacles are more crowded. Thence the spring chain should start the crossing 
at the largest vacant site clusters (VSC). The typical size of such VSC must 
now be determined. Actually the mean number of VSC in a row which obstacle 
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Fig. 3. (color online) Schematic representation of the model used for a quasi-straight spring chain tightly 
bound to a single lattice row. The small full circles represent the lattice sites, the large open circles represent 
the obstacles and the triangles arc for the spring chain nodes. The average number riv of vacant sites involved 
in the largest vacant site cluster is determined through Eq. [6] and Eq. [5] The average spacing m between 
the obstacles on both sides of the largest vacant site cluster is fixed by Eq. [7] 



The normalized probability 
- Cm)" while the probability 



density is fixed to Cm is Lysc = {cmLy — 1) ~ Nm 
to find a VSC with exactly n vacant sites is Cm(l 
for a VSC which size is inferior to n is [1 — (1 — Cm)"]- The probability to find 
a VSC of size n and {Lysc ~ 1) VSC with size inferior to n is proportional to: 



7(n) = LvsC [Cm(l - Cm)"] [1 - (1 - Cm,)"] 



(5) 



The mean size of the largest VSC in the denser row is thus: 



n. 



(6) 



Such a maximum VSC is surrounded by other VSC's that mean size is given 
by: [En<n, '^Cm(l-Cm)"]/[En<n, Cm(l-Cm)"] which for Convenience is denoted 
as (m — 1) with: 



1 (1 - Cm)"" 

m = riy z : — 

c 1 — fl — C 



(7) 



To compute the external force associated with the strongest binding row we 
consider the segment of vacant sites as embedded into a regular lattice 
of obstacles spaced by a mean distance m. Such a mean-field construction is 
illustrated within Fig. [3l where spring chain's nodes (triangles) are bound to 
the lattice sites occupied by the obstacles (large open circles). The array of 
obstacles is assumed to be centro-symmetric, so we ascribe the label to the 
center of symmetry which corresponds to the top of the bulge. We also define 
a new variable n = (1 + n^)/2 for convenience of the notations. Under the 
external applied force r, The force balance sheet, for say the left hand side of 
the chain leads to the set of equations: 
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Fo = -T - 2(xi - xo) 

Fi = -T - {x2 + a;o - 2xi) 

F2 = -T - (X3 + Xi - 2X2) 



n+m) 



Fn-l = -r - (X„ + Xn-2 " 2x„_i) 

Fn = -T - f{Xn) - {Xn+1 + Xn-l - 2x„) 
= -r - {Xn+2 + Xn- 2Xn+l) 

Pn+m—1 T~ {Xn+m ~\~ •^n+m— 2 2Xn+m—l) 

Fn+m fiXn+m) i^n+m+l ~l~ -^n+m— 1 2"^' 

-fri+m+1 (•^n+m+2 ~l~ X^+m 2x^_|_m+l)) (8) 

and in principle the series of equations repeats up to the chain boundaries 
with increment of subscripts. We assume that the mechanical equilibrium is 
satisfied for all nodes j situated in between obstacles. Then we have Fj = 
but for j G [n, n + m, ..n + pm]. For the segment j G [0, n], it is easy to show 
by recurrence that : xj — xq = —Tp/2. Thence the chain shape is parabolic 
between n and —n. For j G [n, n+m], we proceed the same and find x„ = 

—rj{n + j/2) — [Fn + f{xn)]j which fixes the segment end to 

Xn+m = Xn- Tm{n + m/2) - m[Fn + f{Xn)]- (9) 
The same can be iterated once again for j G [n + m, n + 2m] which leads to 

Xn+2m - Xn+m = -Tm{n + 3m/2) - (F„ + Fn+m + f{Xn) + f{Xn+m))m. The 

set of equation on the positions Xn+jm is then: 



Xn+m = Xn - Tm{n + m/2) - f\xn)m 
Xn+2m = Xn+m ~ Tm{n + 3m/2) - [F„ + 

Pn+m + f{Xn) + f{Xn+m)]m 



Xn+pm — Xn+{p-l)m ~ Tm{n + {2p — 1)"^) 

Subtracting the two latest equations yields: 

P ^mXn+pm n/ \ (11\ 

^ n+pm Tm J [Xn+pm ) \^^) 

m 

where ApXn+pm = (a;„,+(p+i)m + Xn+{p-.i)m - "^Xn+pm) IS the discrete Laplacian 
applied to the p subscript. When the entire chain is at mechanical equilibrium 
Fn+pm = for all p. Far enough from the n^-VSC (i.e., the VSC with riy vacant 
sites), the solution for Xn+pm tends asymptotically to a constant Xoo such as 



- m^[Fn+jm + f{Xn+jm)]i^^) 
7=0 
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rm = — /(xoo) and therefore: 

2 arccos (-rm//M) , 47r 
3^00 = ^ cos( + —). (12) 

We can expend linearly Eq.[TT]for the far enough sites such as the displacement 
Xn+pm writes Xn+pm = x^ + Cp and f{xn+pm) = f{xoo) + f'{xoo)ep. Then, at the 
equihbrium Eq. [TT]yields [Apep = —f'{xoo)'mep] and thence ep is an exponential 
function: [ep = eoexp(— ap)] which the exponent a verifies 



a = ±2 ash(y'-/'(xoo)m/2). (13) 

Since the chain displacement is bounded, we are solely concerned with so- 
lutions such as (ap) > 0. The sum of the whole set of equations in Eq. [TOl 
provides another relation between r and the nodes position 

Xn+pm.} the 

condition that Fn+pm = for all p: 

p-i 

Xn - Xn+pm = m[Tp{n + pm/2) + "^{p ~ j)f{Xn+jm)], (14) 

i=o 

which after expanding f{xn+jm) as a Taylor series around Xoo and keeping 
only the terms linear in p provides us with an equation which relates r to eo: 

"^o r f'{Xoo) , ^of"{xoo) I 4f"'{xoo) , , 



(n-m/2)'(l-e-") 2(1 - e'^") 6(1-6" 



The critical chain configuration is reached when the Hessian associated with 
Eq. [11] has a singular eigenvalue. This allows us to determine the critical 
value for eo. Actually we found that finding the Hessian singular eigenvalue 
is equivalent to find the maximum of Eq. [15] for r with respect to eo. The 
solution for the critical bulge is then: 

^ / (3^00) \j f '{Xoo)"^ 4/"'(Xoo)/'(3^oo) 3(]^„e-3a^(]^lg-~ 

'° = "^^^ " '"'"^ 2r(xoo)(l-e-2") • 

(16) 

Combining the solutions for Eqs. [12], [131 IE] and [16] allows us to determine 
the maximum pinning force associated with N^, the number of obstacle in 
the denser row. For this reason, we denote such a maximum as r(iVm). The 
set of equations giving r(A^m) can be solved recursively. Starting with a small 
enough trial solution for r = ro, we compute the corresponding quantities Xoo 
and a from Eq. [T2] and from Eq. [T3] Then eg is derived from Eq. [TH] and the 
corresponding value of r from Eq. [15] If the so obtained quantity is larger than 
the initial value Tq then the latter is incremented and we proceed the same up 
to find identical values for r and To. The end result gives the required T{Nm) 
to a precision fixed by the trial solution increment. 
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Fig. 4. (color online) External force required to drag a spring chain of length Ly over a distance Lx, for 
a pinning strength /j;/ = 0.005 and a cutoff w = 0.1 in (a) and w = 0.5 in (b). The symbols represent 
the data obtained through the simulations described in Sec. [21 for different lattices (see figures legend). 
The continuous lines correspond to the predictions made through the theory detailed in Sec. \3\ for same 
parameters as those used in simulations. Colors of symbols and lines correspond one to one. 
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Fig. 5. (color online) External force required to drag a spring chain of length Ly = 1000 along a distance 
Lx = 100, for different obstacle pinning strengths (see legend) and the same interaction range w = 0.5. 
The different symbols represent the simulations data and the continuous lines correspond to the prediction 
made through the analytical theory detailed in Sec. |3] for same parameters. 

The critical pinning force Tc of the random lattice is approximated by averaging 
T{Nm) over A^^: 

n = Y.P{NmMN^), (17) 

where (3{Nm) has been given in Eq. HI The previous theory is compared to sim- 
ulations data in Figs. H] (a) and (b) and in Fig.[5]for different lattice dimensions, 
different pinning forces, varying w and J'm- A quantitative agreement has been 
obtained between theory and simulations, although no adjustable parameters 
are involved. According to same type of comparisons but for larger /m values, 
the previous analytical work proves relevant for Jm smaller than roughly 0.03. 




10 




position of the chain / lattice X-axis 

Fig. 6. (color online) Strongest pinning configuration for fj^j = 0.03, Ly = fOOO, = fOO, a density 
Cs = 9 % and an interaction cutoff w = 0.5. The inset shows a chain segment the nodes of which are marlced 
with triangles. Only obstacles close from the string have been reported as circles. 

The theory predictions worsen for cases where the critical configuration crosses 
few lattice rows. As an example, for = 0.03 in Fig. [5l the VSC theory is 
found to become less accurate for low densities. In such conditions, the string 
profile at the depinning transition corresponds to a kinked shape, as shown 
in Fig. [6l different from the wavy profile shown in Fig. [1] and from the quasi- 
straight ones shown in Figs. [2] (a) and (b). From Fig. it can be seen that the 
discrepancy increases as the density decreases while the theoretical predictions 
remain accurate for more concentrated obstacle distributions. As /m increases 
above 0.03, the deviation between theoretical predictions and simulations data 
is shifted toward higher densities. In this range of parameters, the system 
undergoes a bifurcation, not treated in the present work. 

Other authors [?,?] noticed that the maximum pinning force of a random lat- 
tice was dependant of the drag distance. Concerning the weak pinning points 
studied here, it is thus of some interest to explore the variation of the critical 
drag force with lattice dimensions. The critical drag force was found to vary 
proportionally to [ln(L2;)]""= where varies with all parameters. For instance 
when w = 0.5, /a/ = 0.01 and Ly = 1000, we found = [0.45 - 0.036 ln(c,)]. 
One can thus conclude that the dependence is very weak since a frac- 
tional power of a logarithm is a rather wise function. The maximum drag 
force depends not only on but also on Ly. The critical drag force varies 
proportionally to a constant plus the function [ln(Ly)]~"!' where ay decreases 
when Cs increases. For instance when w = 0.5, /m = 0.01 and L^ = 4000, we 
found ay = [0.26 — 0.631n(cs)]. According to direct simulations, the Tc depen- 
dence in Ly seems to weaken when J'm becomes large enough for the critical 
configuration to differ from the quasi-straight line as in Fig. [6] or Fig. [H 

For fixed lattice dimensions, the adjustment of a density power law for the 
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theoretical critical drag force allows to establish some comparison with the 
SSH analytical theory. The density power law fit is found very close from a 
linear variation and may even overpass slightly the unitary exponent in some 
situations depending on the lattice geometry. For instance, from Fig. |l](b), we 
worked out by a curve fitting with the form Tc oc (c^), an exponent t] = 1.13 
for Ly = 4000 and = 25 whereas for Ly = 500 and = 1000, t] = 0.946 
was obtained. The effective density exponent rj is therefore size dependent. 
A density exponent close from unity corresponds to the FMS theory for the 
hardening of concentrated solid solutions. The CRSS linear dependence in the 
solute concentration was also noted in atomistic calculations on the model 
solid solution Ni(Al) [?]. 

3.2 Extension to broader cutoff' 

To extend the theory developed previously for a short cutoff w, we first con- 
sider the case w ~ 1, which allows us to treat the spring chain interaction 
with only two rows: the row at the back and the next nearest one. Such a 
critical configuration can be seen from direct simulations as reported in Fig. 
H] (b). To compute Tc, we still employ the same model as depicted in Fig. [31 
However, the interaction with the next nearest row may weaken the pinning 
on the back row. One ascribes to each site j of the next nearest row, a force 
g{xj) = Cvf{xj — 1) where c„ is the density of obstacles in the next nearest row. 
For the situation shown in Fig. [31 the force balance sheet writes as follows: 

Fo = -T - g{xo) - 2(xi - a;o) 

Fi = -T- g{xi) - {x2 + xo - 2xi) 

F2 = -T - g{x2) - (x3 + xi- 2x2) 

Fn-l = -r - g{Xn^l) - {Xn + Xn-2 " 2x„_i) 

Fn = -T - g{Xn) - f{Xn) " {Xn+1 + X^-l - 2x„) 
Fn+1 = -T - g{Xn+l) - {Xn+2 + X^ - 2x„+i) 

Fn-i-m—1 fi'('^n+m— 1) (■^n+ni ~l~ ■^n+m— 2 2Xn+m—l) 

Fn+m T~ g{Xn+m) f^Xn+m) {Xn+m+1 ~l~ Xn+m~l 2Xn+m) 

-fn+m+1 T fi'('^n+m+l) (•^n+m+2 ~l~ -^n+m 2Xn+m+l) 

(18) 

To estimate the force exerted by the next nearest row upon the spring chain, 
the function g{xj) is approximated with a step function: g{xj) = g{xo) if 
j E [0; n - (m - l)/2[ and g{xj) = g{xn+pm) if j e [n + pm - {m - l)/2; n + 
pm + (m+ l)/2[. We then assume that the mechanical equilibrium is achieved 
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for all nodes excepted those aligned with some obstacles of the back row. This 
leads then to the following set of equations: 



T = -g{xo) - 2(xi - Xq) 

T = -g{xo) - (x2 + a;o - 2xi) 

r = -g{Xn) - {Xn + Xn-2 - 2x„_i) 

T f{-^n) (-^n+l ~^ -^n— 1 Sx^) 

r = -g{Xn) - {Xn+2 +Xn- 2x„+i) 



n+m— 2 



- 2x 



n+m—l 



-^n+m fiS^n+m) diS^n+m) ~l~ •^n+m—l '^•^n+rri) 

T = —g{Xn+m) ~ {Xn+m+2 + Xn+m ~ 2x„+m+l), (19) 

and the equations repeat up to the chain boundaries by incrementing sub- 
scripts. By applying recurrence, it is possible to reduce the previous set of 
equations to a smaller one, concerning only the regular array of obstacles in 
the back row: 



p-i 

Xn+pm = Xn+{p-i)m " rm{n + {2p - l)m/2) - m\Y^\Fn+jm + 



p-2 

+m'^g{xn+jm) + g{xo){n - m/2) 

j=0 

-m giXn+{p-l)m)- 



— 1 



g (^^n+pm) 



(20) 



Subtracting the equation for rank (p — 1) from the one at rank (p) gives: 



n+pm ■ 



^m-^n+pm TTfl n^'fi^^n+pm) 



m 



9{Xn+{p+l)m) + 9{Xn+{p-l)m) 



3m2 + 1 



9i.Xn+pm), (21) 



while for p = 0: 



Fn = Xn - Xn+m " Tm{n + m/2) - m[f{xn) + 9{xo){n - m/2)] 



rn? — 1 



8 



-g[Xn+ni) - m- 



3m — 1 



-9[Xn) 



(22) 



When the equilibrium is achieved Fn+pm = for all p, which allows us to 
;he asymptotical solution of E 

verifies: 



deduce the asymptotical solution of Eq. [2T] as a constant Xn+pm which 



rm = -/(xoo) - mg{xoo)- 



(23) 
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Fig. 7. (color online) Same as in Fig|3]but for a cutoff w = 1 and a pinning strength = 0.001. 

The solution to Eq. [23] is the positive real root of a third order polynomial 
equation: 



1 + mcv 
with Aa 

and Br^, 



cos arccos 



-5. 



27 X , 4:^T^ 



— 1 — mCt, — 



1 



+ 



1 + mCy 
3W^" 



^2(1 + mc^ 



w(l + mCy) 



(24) 



If (1 — Xoo) > w then we can set c„ = in the previous equation which 
leads to Eq. [121 valid for small w. We now expend linearly Eq. [21] around Xoo 
to determine an approximation of nodes position as Xn+pm = Xoo + e^. The 
equation on is: 



-m/'(xoo)ep 



3m^ 



— 1 



4 



g'{xoo) ep+i + ep_i . (25) 



The solution is an exponential function ep = eoexp(— ap) with a dispersion 
relation: 



/ 2 - mfixoo) - (3m^ + l)g'{xoo)/A ^ 
[ 2 + ^.9'(Xoo) ^' 



a = ±ach( - — Y^v'^-^^- ). (26) 

At this stage, it is of some interest to work out the maximum of the spring 
chain position. For the segment situated along the larger hole (see Fig. [3]), 
j G [0,n — (m — l)/2[, the set of very first equations in Eqs. [19] leads to: 



rn 



rrr — 1 



Xfi Xq 



-9[Xn) 



\n — m' 



74 + 1/4]. 



(27) 
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Then xq can be expressed as a function of x„ since xq is actually the positive 
root of a third order polynomial: 



Xo = l + 2w]^ 
with Aq 



I Aq + w 



cos 



arccos 



Co 



27AI 



2Ao\ iAo + 



w 



/3 + 47r/3 



n 



mV4 + 1/4 



and Co = Xn — I 



Tfi m 
^ + — 



'-9{Xn)- 



(28) 



We also need to express the first derivative of xo against Xn which according 
to Eq. EH gives: 



dxo 



1 + r^g'jxn) 

[^2 _ (^2 _ l)/4]c,'(a;o)/2' 



(29) 



The sum of the equations in Eq. [20] from rank 1 to rank p, taken sufficiently 
large, leads to an equation which the linear term in p is: 



r 



-9{xo) + 



[n — m/2] 



5m 



eo- 



(/'(a^oo) +mg'{x^)) 



(/"(Xoo) + mg'\x^)) 3 (/"'(Xoo) + mg"\x^)) 



2(1 



-2a^ 



6(1 



(30) 



The latter equation is similar to Eq. [15] obtained for the tightly bound chain, 
but includes the interaction with the next nearest row. The maximum r in 
Eq. [30] against eo corresponds to the critical strength which provides a tran- 
scendental equation on eo: 



[n-m/2\g'{xo) 



dxo 



if'ixoo) + mg'{xoo)) (/"(a^oo) + rng'^Xoo)) 
eo — 



1 - e-" 
,{f"{x^)+mg"'{x^)) 



l-e 



-e, 







2(1 



-3a\ 



(31) 



where Xn = x^o + Co; and dxo/dx^ are given in Eq. [23 and Eq. [23 This 
complete our computation of the maximum pinning force associated with 
and Cy. The corresponding value of r is therefore related to A''^ through Eq. [6] 
as well as to the number of obstacles iV„ = CvLy. The maximum force is now 
considered as a function of both Nm and A''^ and it is denoted by T{Nm, N^)- 
The probability to find a couple of rows which actually consists of a back row 
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Fig. 8. (color online) External force required to drag a spring chain of length Ly = 200 over a distance 
Lx = 1500, for different obstacle pinning strengths (see legend) and the same interaction range w = 1. 
The different symbols represent the simulations data and the continuous lines correspond to the predictions 
established through the analytical theory detailed in the text. 

with Nm obstacles and a front row with A'^ obstacles as being the strongest 
configuration among rows is written as: 

^(iV^,iV,) = [p{Njp{N,)] [1 - E p(iVi)p(iV2)]^^-\ (32) 

r{NuN2)>TiN^,N.,) 

where the function p{N) is the binomial given in Eq. [31 Thence the average 
critical depinning is given by: 

Tc= r{N„^,N,)e{N„,,N,). (33) 

The same statistical treatment can be applied to the case w < 0.5 and yields 
the same results as presented in Sec. 13.11 since the critical force given in Eq. 
[15] is then independent of A''^,. 

In Fig. [7] and Fig. [HI the theoretical predictions are compared to simulations 
data. It is worth noticing that the critical force variation against density cor- 
responds to a more convex curve than for shorter cutoff w < 0.5. An analysis 
in term of power law fit leads to an exponent smaller than unity. For instance 
in Fig. [7], for Ly = 500 and = 1000, the effective density exponent is around 
a = 0.65. The effective density exponent of Tc is found to decrease when w 
increases (see further Sec. 13. 3p . In Fig. [7], the computations have been per- 
formed for different aspect ratios of the lattice and show the same trend as 
in Fig. [H for a shorter w, the critical force density exponent increases with 
Lj. and decreases with Ly. Although a satisfactory agreement is obtained for 
the different lattices, we notice that the pinning force Jm is smaller than in 
Fig. [H Actually the field of validity for the VSC theory is narrower for w = 1 
than for w < 0.5. The limit of application for the theory decreases as the 
interaction cutoff increases. Such a limit also varies with the lattice dimen- 
sions: it decreases when Ly increases and when L^ decreases. For instance. 
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with Ly = 200 and = 1500, the theory proves to be efficient as seen from 
Fig. [8] where the comparison has been performed for different pinning forces 
the maximum of which is J'm = 0.006. The same computations carried out for 
Ly = 2000 and = 500, keeping constant both w and /m yield much worse 
results, in particular for the low densities. This can be understood comparing 
the critical proffies for both geometries. Whereas the chain profile is quasi- 
straight for the former the critical profile is wavy for the latter, i.e., similar to 
the one seen in Fig. [TJ As for short cutoffs, the system undergoes a bifurcation 
passing from a quasi-straight critical proffie to a wandering one. The change 
in critical proffie occurs for smaller J'm with w = 1 than for w < 0.5. 



3.3 Extension to intermediary and still broader cutoffs 



In previous computations, we developed the pinning force function as a Taylor 
series around Xoo, i.e., the asymptotic solution for the chain nodes position. 
For intermediary cutoff, namely 0.5 < w < 1, we must be aware that such a 
development cannot be used to approach the non-analytic force function since 
the cutoff occurs right in between two rows. Then Xoo is assumed to remain far 
from the cutoff interaction with the next nearest row. The asymptotic position 
Xoo is thus given by Eq. [12] and a by Eq. [131 The top of the bulge is assumed 
to be situated above the cutoff distance from the next nearest row. Then the 
equation on xq is same as in Sec. 13.21 and can be solved analytically as shown 
in Eq. [28l The equation relating r to cq must be rederived. The sum of the 
equations in Eq. [20] from rank 1 to rank p, taken sufficiently large, leads to an 
equation which the linear order in p is now: 



_ _ / \ ^ \ f (^oo) 2 / (^oo) , 3 / i^oo) 

+ ^9{Xn+jm)\: (34) 

i<ib 

where is such as {xn+it,m > I — w). The quantity lb corresponds to the 
segment of the bulge which overpasses the interaction cutoff with next nearest 
row. Such equation for r holds provided that lb remains small compared to 
Ly/m. 

In Fig. [9] the predictions from the present development are compared to sim- 
ulations data for w = 0.7. We reported the results obtained from simulations 
with different obstacle densities. For some densities, we carried out the simula- 
tions with different random distributions. The simulations data are distributed 
around the theoretical predictions which con&ms the robustness of the theory. 

The previous computations can also be extended to an interaction cutoff larger 
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Fig. 9. (color online) External force required to drag a spring chain of length Ly = 1000 over a distance 
Lx = 200, for a pinning strength /m = 0.001 and different interaction range ui (see legend). The different 
symbols represent the simulations data and the continuous lines correspond to the prediction made through 
the theory detailed in the text. 

than the unit cell parameter (i.e., w > 1). In the theory, the function g{xj) = 
Cvf{xj — 1) is substituted with g{xj) = c^f{xj — 1) + h{xj) where h{xj) = 
Csf{xj — 2) accounts for the pinning force due to the still further row, i.e., 
the second next nearest one. The top of the bulge may interact with such a 
row meanwhile the back of the chain is assumed to remain far from it. Then 
the equations found in Sec. 13.21 for Xoo and a still hold and xq is again given 
by Eq. [27| but substituting the function g{x) with the expression g{xj) = 
Cvf{xj — 1) + h{xj). The equation on r follows: 
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3<h 



where corresponds to the top of the bulge overpassing the distance cutoff 
with the second next nearest row, defined as [xn+i^m > 2 — w). In Fig. [9l 
the theoretical computation for against the obstacle density is compared 
with simulation data for {w = 1.1) and for different random distributions. 
Although the theory could certainly be improved by accounting for the effect 
of the nearest row at the back of the chain, the agreement with simulations 
data proves quantitative. 

The results obtained in both cases, w = 0.7 and w = 1.1 show that Tc decreases 
with w. This was confirmed by other simulations performed for still larger w, 
up to w = 2. Such a decreases is opposite to the one predicted within the 
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MNL theory for SSH. Moreover according to the MNL theory, the density 
exponent of Tc is fixed to t] = 2/3, whereas according to the VSC samphng 
theory such an exponent depends on lattice dimensions (see Sec. 13.21) and w, 
as welL For instance, in Fig. [HI a power law fit on the theoretical predictions 
yields and effective exponent 77 = 0.945 for w = 0.7 while for u; = 1.1 we 
obtained r] = 0.827. The variations 77 with lattice dimensions, along with 
its value larger than 2/3 and its w dependence demonstrate that the MNL 
theory does not apply to the systems studied here. The same conclusion holds 
for other SSH theories predicting a constant density exponent of the critical 
depinning threshold. 



4 Summary and perspectives 

In the present study, we addressed a discrete version of a paradigmatic prob- 
lem, namely the depinning threshold of an elastic string on a random substrate. 
For a planar distribution of weak pinning points, a theory was devised to com- 
pute the applied force required to drag the one-dimensional elastic manyfold 
over a disordered potential landscape with various aspect ratios. The the- 
oretical predictions were found accurate provided the critical configuration 
remains close from a quasi-straight line bounded between two lattice rows. 
The strongest pinning configuration was shown to hinge on the denser lattice 
rows in which the largest vacant site clusters (VSC) are bounded in size. Such 
maximum VSC correspond to the weakest lattice defects on which the critical 
depinning proceeds. The mean size of the critical VSC is determined through 
an expression involving only lattice dimensions and the overall planar obstacle 
density Cg. The theory allowed us to account for the finite lattice size effects 
in a quantitative manner. The typical variations of the critical applied force 
against the chain length and the drag distance both yield a logarithmic power 
law. For a fixed lattice geometry, and an interaction cutoff inferior to half the 
lattice cell parameter the pinning strength was found close to being linear in 
Cg. The effective density exponent rj of the depinning threshold was found to 
depend on the lattice geometry and to decrease as the interaction cutoff w 
increases. 

In some atomistic studies bearing on dislocation in a model Ni(Al) solid so- 
lution, i.e., with a rather marked size effect [?], the maximum force exerted 
by isolated solute atoms on a dislocation segment was found of the order of 
0.04 nano-Newton (nN). The elastic contribution to the line tension of a screw 
dislocation, computed within the isotropic elastic theory [?] is F = rsln{R/b) 
where b is the Burgers vector, R is the inter-dislocation spacing and where 
the pre- logarithmic factor is given by F^ = /ifc^(l + z/)/47r(l — z/), with /i as 
the shear modulus and u as the Poisson coefficient. The distance R is related 
to the dislocation density pd hj R = l/^/pd which for a standard density in 
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deformed metals pd = 10^^ m~^, gives R = 1 fim. Then, with fiNi = 74600 
MPa and u = 0.28, the screw dislocation line tension in Ni is F = 6.1 nN 
which gives a ratio between the obstacle pinning strength and the line tension 
of Im = 0.007, of the order of /m studied in the present work. A solute atom 
dilation smaller than the one for Al in a Ni matrix or a smaller dislocation 
density could even yield smaller Jm- In regard of the simplicity of the elastic 
manyfold model, the present study requires to be extended to a more realistic 
model involving some of the atomistic details as for instance: (i) a position 
dependant obstacle strength, (ii) a dissociated core, and (iii) mixing of differ- 
ent types of obstacles. Such extensions would allow to address longstanding 
issues in SSH theory. 
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